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ABSTRACT 

In low collisionality plasmas heat flows almost exclusively along magnetic field lines, and the con- 
dition for stability to convection is modified from the standard Schwarzschild criterion. We present 
local two and three-dimensional simulations of a new heat flux driven buoyancy instability (the HBI) 
that occurs when the temperature in a plasma decreases in the direction of gravity. We find that 
the HBI drives a convective dynamo that amplifies an initially weak magnetic field by a factor of 
^20. In simulations that begin with the magnetic field aligned with the temperature gradient, the 
HBI saturates by rearranging the magnetic field lines to be almost purely perpendicular to the initial 
temperature gradient. This magnetic field reorientation results in a net heat flux through the plasma 
that is less than 1% of the field-free (Spitzer) value. We show that the HBI is likely to be present 
in the cool cores of clusters of galaxies between ~ 0.1 — 100 kpc, where the temperature increases 
outwards. The saturated state of the HBI suggests that inward thermal conduction from large radii 
in clusters is unlikely to solve the cooling flow problem. Finally, we also suggest that the HBI may 
contribute to suppressing conduction across cold fronts in galaxy clusters. 

Subject headings: convection — instabilities — MHD — galaxies: clusters — plasmas 



1. INTRODUCTION 

For thermally stratified fluids, stability to convection 
is guaranteed if the entropy decreases in the direction 
of gravity. In dilute astrophysical plasmas, however, the 
mean free path is large compared to the gyroradius, and 
heat condu ction is anisotrop ic with respect to the mag- 
netic field (Braginskii 1965). In the low collisionality 
regime, the convective instability condition is modified, 
and plasmas in which the temperature inc reases in the d i- 
rection of gravity are buoyantly unstable (jBalbus Il2000f ) . 
This instability, called the magnetothermal instability 
( MTI), has been simulated in tw o and three dimensions 
bv lParrish fc Stone I 12001 . 

In this Letter, we study a related in stability of low col- 
lisionality plasmas. iQuataertl (|2007l ) has shown that, in 
the presence of a background heat flux, anisotropic con- 
duction drives a buoyancy instability when the temper- 
ature decreases in the direction of gravity. The growth 
rates of this heat flux driven buoyancy instability (HBI) 
are verified using two-dimensional (2D) simulations in 
fj2] In <j3] we determine the nonlinear saturation of the 
HBI in three dimensions (3D) and show that it drives a 
magnetic dynamo while greatly suppressing the net heat 
flux through the plasma. Finally, in $4] we apply the HBI 
to the intracluster medium and discuss its implications 
for cooling flows and cold fronts. 

2. 2D SIMULATIONS 

2.1. Equations and Method 

We solve the usual equations of magnetohydrodynam- 
ics (MHD) with the addition of anisotropic thermal con- 
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duction. The MHD equations in conservative form are 
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where the symbols have their usual meaning. The total 
energy E is given by 
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where e = p/("f — 1). Throughout this paper, we assume 
7 = 5/3. The anisotropic electron heat flux is given by 

Q = -xebb ■ VT, (6) 

where xc is the Spitzer conductivity (jSpitzer Ill962f ) and 
b is a unit vector in the direction of the magnetic field. 

We follow the HBI into the nonl inear regime using 
the c onservative MHD code Athena (Gardine r fc Stone I 
I2005D with the addition of anisotropic conduction along 
magne tic field lines as described in iParrish fc Stone I 
(2005). This methodology has been extensively utilized 
and tested previously in the study of the MTI. In all 
simulations, g = 1. The 2D simulations use a linear 
temperature gradient with fixed temperature boundary 
conditions at the top and bottom of the domain. 
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Fig. 1. — Snapshots of the magnetic field in a 2D simulation 



initialized with a single mode perturbation having k x 



The 



HBI drives the initially vertical field to become largely horizontal. 
The units of time in Figs. 1-4 are such that the dynamical time 
(gdlnT/dz)' 1 / 2 » 1.4. 

2.2. Physics of the HBI 

The physical origin of the HBI can best be under- 
stood by examining magnetic field snapshots from a 2D 
simulation, as shown in Figure Q] (see also Fig. 1 of 
iQuataert Il2007fh This 2D simulation was run at a res- 
olution of (100) 2 . The initial temperature profile was 
T(z) — Tq(1 + z/2), the temperature scale-height at the 
midplane is 2, and the size of the domain is (0.1) 2 ; note 
that the simulations arc local. The pressure and den- 
sity were chosen so that the atmosphere was in hydro- 
static equilibrium with p ~ p ~ 1. The magnetic field 
is chosen to be weak initially and purely vertical with 
Bo/ '{Ait) 1 / 2 = 5 x 10~ 4 , so that magnetic tension forces 
are negligible. The purely anisotropic thermal diffusiv- 
ity is k = xc'T/P = 10~ 2 , where k has units of a dif- 
fusion coefficient, i.e., cm 2 s _1 . With these parameters, 
the conduction time for small-scale perturbations is much 
less than the dynamical time, which is the limit of fastest 
growth for the HBI (and MTI). 

Figure [1] shows that perturbations with non-zero k x 
and k y generate converging and diverging field lines. The 
heat flux follows these field lines leading to (conductive) 
heating and cooling of the plasma. In a plasma with 
dT/dz > 0, a fluid element displaced upwards is thus 
heated by the background heat flux, causing it to rise 
further and become buoyant. Note that the magnetic 
field snapshots in Figure [1] look very similar to snapshots 
for the MTI rotated by 90°. 

For a weak vertical magnetic field, the growth rate of 
the HBI in the limit of rapid conduction is 
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For the parameters of our simulation, we predict a growth 
rate of 0.5 for k x = k y . Measurement of the growth 
rate in our single mode simulation verifies this prediction 
to within 1%. In the 2D simulation, we find that the 
magnetic energy is amplified by a factor of ~ 52 during 



the course of the run. We defer a detailed discussion of 
the nonlinear saturation to the 3D simulation. 

3. NONLINEAR SATURATION IN 3D 

Only in 3D can we accurately explore the saturation 
of the HBI since the nature of convection differs signifi- 
cantly in 2D and 3D and th e anti-dynamo theorem ap- 
plies in 2D (|Cowling | [l934). Using 3D simulations, we 
now quantify the amplification of the magnetic field by 
the HBI and the resulting heat flux through the plasma. 
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Fig. 2. — Time evolution of the volume-averaged magnetic energy 
density and kinetic energy density in the HBI-unstable region of 
run Rl. The magnetic energy is amplified by a factor of ~ 300 and 
is in rough equipartition with the subsonic convective motions. 

The computational set-up we use is nearly identical to 
the st ratified box presented in §2.3 of iParrish fc Stonel 
([2001 : we refer the reader there for details. The box 
consists of three vertical layers of equal size in which 
the central region has a linear temperature profile with 
dT/dz > and pure anisotropic conduction along mag- 
netic field lines, such that it is unstable to the HBI. The 
top and bottom layers are isothermal atmospheres (expo- 
nential pressure and density profiles) and act as a buffer 
to the penetrative convection that takes place. These 
regions have isotropic conductivity and are stable to the 
HBI. For all of the simulations in this paper, the isotropic 
conductivity in the top and bottom layers is equal to the 
parallel anisotropic conductivity in the middle region. It 
is important to note that our box size is again small 
compared to the temperature scale height (0.2 and 2, re- 
spectively), so that our simulations are local in nature. 
Table [T] gives the initial magnetic field and conductiv- 
ity for our 3D runs; we will primarily focus on Rl, our 
fiducial simulation. 

Figure [2] shows the amplification of magnetic energy as 
a function of time for Rl. The net amplification can be 
defined by 

A(B 2 ) s i*3*L, ( 8 ) 

(-D 2 )init 
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TABLE 1 
Table of Nonlinear Runs 



Run 


Bo/(47r) 1 /2 


re II 


A{B 2 ) 


RMS Mach 


(8 b) 


/ = Q/Q 


Rl 


icr 5 


5 x 10" 3 


327 


5.1 x 10" 4 


5.8° 


0.51% 


R2 


icr 4 


5 x icr 3 


36.4 


7.7 x icr 4 


17.2° 


7.7% 


R3 


icr 5 


1.5 x 10" 2 


384 


5.1 x 10" 4 


3.8° 


0.31% 



which for the fiducial case is A(J5 2 ) ~ 327. In addition 
to amplifying the field, the HBI drives vertical convec- 
tive motions. These remain subsonic with an RMS Mach 
number of 5.1 x 10~ 4 ; the kinetic energy is ~ 10 times 
the magnetic energy at the end of Rl. 

The evolution of the average angle of the magnetic field 
with respect to the horizontal in the unstable layer is 
particularly interesting; this is defined by 

« s ( si °"'(w))- (9) 

Figure [3] shows (6 b ) for run Rl. The buoyant motions 




Fig. 3. — Evolution of the average angle of the magnetic field with 
respect to the horizontal for the fiducial simulation Rl. The HBI 
drives the initially vertical field to become almost purely horizontal 
with an average angle of only 5.8°. 

of the HBI saturate by driving the magnetic field to 
be largely perpendicular to the background temperature 
gradient. This evolution is in stark contrast to the MTI, 
which preferentially aligns the magnetic field with the 
background temperature gradient. For Rl, the magnetic 
field evolves from an initially vertical field (90°) to a field 
with an average angle of (6b) ~ 5.8°. As one might ex- 
pect, this dramatically decreases the heat flux through 
the plasma, as is shown explicitly in Figured] It is use- 
ful to compare the heat flux in the saturated state of the 
HBI with the heat flux that is present in the initial state, 

Q = -*\\(^) • (io) 
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horizontal (Fig. [3]), 
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Fig. 4. — Time evolution of the horizontally-averaged vertical 
heat flux at the midplane in the fiducial simulation Rl. The total 
heat flux (solid line) is divided into the dominant conductive flux 
(dashed line) and the advective flux (dotted line). Note that the 
total heat flux and conductive flux are essentially collinear. These 
fluxes are normalized to the fiducial heat flux, Q, that would be 
present for vertical magnetic field lines (eq. 1101 ). 

layer plummets. We define the conduction efficiency by 
/ = Q/Q- The saturated state of the HBI has a very low 
heat flux of / = 5.0 x 10 -3 for Rl; most of this flux is car- 
ried by conduction, with only a small contribution from 
convective motions (Fig. 01 . We will discuss the obvious 
implication of this result for cooling flows shortly. 

For a weak initial magnetic field, the saturation of the 
HBI is primarily due to th e reduction in th e mean angle 
of the magnetic field (6 b)- iQuataertl (|2007t ) showed that 
for nearly horizontal fields the maximum growth rate of 
the HBI is given by 



M ~ (gdlnT/d*) 1 / 2 ^} 



(11) 



At the conclusion of run Rl, the growth rate has de- 
creased by a factor of ~ 10, to co ~ 0.07, significantly 
reducing the efficacy of the instability. Meanwhile the 
growth of the magnetic field increases the effects of mag- 
netic tension. The Alfven frequency (a lengthscale de- 
pendent quantity) ranges from 0.01 to 0.7 from the size of 
the domain down to the grid scale, respectively. Because 
the Alfven frequency is comparable to the dynamical fre- 
quency on small-scales, and to the growth rate on all but 
the largest scales, tension will further inhibit the growth 
of the HBI. Although we cannot rule out modest addi- 
tional decrease in (6b) and the heat flux on timescales 
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longer than we have simulated, such evolution will take 
progressively longer as the growth rate decreases further 
(eq. and is unlikely to modify the magnetic and 

kinetic energy, which have saturated (Fig. [2j. 

For completeness we ran several additional simulations. 
Run R3 has a conductivity three times higher than the 
fiducial value. However, since Rf is already in the limit in 
which the conduction time for small-scale perturbations 
is rapid compared to the dynamical time, little change is 
seen from the fiducial simulation. Run R2 has an initial 
magnetic field one order of magnitude larger than the 
fiducial run. Due to magnetic tension effects on the scale 
of this local simulation, R2 saturates earlier and with less 
magnetic field amplification. The magnetic field is still 
primarily horizontal upon saturation and the heat flux is 
strongly suppressed, although by less than in the fiducial 
simulation. In a global simulation, the saturated state of 
R2 would be unstable on large-scales, which would lead 
to additional field amplification and rearrangement. 

4. APPLICATIONS AND DISCUSSION 

The HBI is predicted to be present from ~ 0.1 — 100 kpc 
in the intracluster medium of galaxy clusters, where the 
observed temperature increases outwards. For concrete- 
ness, consider the cluster AI795 which has an estimated 
viria l mass of 1.2x 10 15 and a scale radius of r s = 460 
kpc (jEttori et al.ll2002HZakamska fe Nara van 2003). The 
temperature increases from ~ 2 keV at 10 kpc to ~ 7.5 
keV at 100 kpc. Given these parameters, the growth rate 
of the HBI at R — 50 kpc in the limit of weak fields is ~ 
3x 10~ 16 s _1 , i.e., a growth time of ~ 10 s yr. This should 
be compared to the wavelength(A)-dependent Alfven fre- 
quency of ~ 5 x 10~ 17 (X/R)- 1 (B/lfiG) s" 1 given the 
measured electron density of n e ~ 0.02 cm~ 3 at 50 kpc 
and an assumed B ~ l^iG. Thus tension has a small ef- 
fect for A > O.li?. Because the electron mean free path is 
~ 0.05i?, the thermal conductivity is sufficiently high for 
perturbations to grow on the dynamical time (eq. [7]). 

Our simulations presented in §3 were local (L z < 
H), but their nonlinear outcomes provide a clear les- 
son for galaxy clusters. So long as magnetic tension 
is not dominant, the HBI rapidly and efficiently re- 
orients the magnetic field lines to be perpendicular to 
the background heat flux, resulting in a net heat flux 
through the plasma that is much less than the field- 
free Spitzer value. Our r esults are significantly differen t 
from the prediction s of | Naravan fc Medyedev I (2001), 
IChandran fc Maronl (|2004h . and iLazarianl (|2006h . all of 
whom highlight the effects of field line wandering on 
the effective thermal conductivity of a plasma. In - 
stead, however, o ur results, and those of lBalbusl (|2000f) . 
iParrish fc Stone I (|2007l ). & iQuataert I (j2007j), demon- 



strate that the dynamical coupling between anisotropic 
thermal conduction, the heat flux, and the magnetic 
field geometry is crucial. For HBI-unstable plasmas, the 
net result is that heat is transported far less effectively 
than previously assumed. These conclusions make the 
cooling flow problem in clusters more severe since the 
HBI strongly suppresses the transport of energy from the 
outer parts of clusters (r > 100 kpc) to the cooling core. 
An alternate mechanism is thus necessary to heat the 
cluster plasma, e.g., feedback from a central AGN. Sev- 
eral additional effects should be accounted for before this 
conclusion is considered definitively established. First, 
cosmic rays from a central AGN modify the stability 
prope rties of intracluster plasma ([Chandran fc Dennis! 
2006). The interplay between the cosmic rays, which 
tend to drive the system MTI unstable, and the plasma, 
which by itself is HBI unstable at intermediate radii, may 
be interesting and complex. In addition, turbulence pro- 
duced by galactic outflows and galaxy wakes could mod- 
ify the magnetic field geometry established by the HBI. 
Finally, given that tension restricts the HBI to scales 
> O.li?, global simulations are clearly required. 

An additional potential application of the HBI in clus- 
ters is to the suppression of conduction across cold 
fronts. Th e canonical cold front is found in the clus- 
ter A3667 ()Vikhlinin et alj [2001a. b); our numerical val- 
ues are drawn from this example. The cold front is an 
abrubt change in plasma properties, with high tempera- 
ture, low density, low pressure plasma lying outside low 
temperature, high density, high pressure plasma. The 
temperature gradient is thus opposite to both the pres- 
sure gradient and the cluster's gravity — an HBI unstable 
situation. For the gradient scale length, we take the 5 kpc 
upper limit from Chandra measurements and estimate an 
HBI growth rate of u> s=s 3.8 x 10~ 15 s _1 , corresponding 
to a growth time of 8.5 Myr. For a 1 /iG field, magnetic 
tension is modest for perturbations with wavelengths of 5 
kpc, suggesting that the HBI can indeed grow. Although 
the HBI cannot account for the origin of cold fronts, it 
may help account for their survival. Namely, the HBI 
will drive the magnetic field to be horizontal in the cold 
front, suppressing heat transport across the cold front 
and thus maintaining the temperature gradient in spite 
of the nominally short conduction times. 
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